Could two NMSSM Higgs bosons be present near 125 GeV? 
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Data from the ATLAS and CMS collaborations [THS] provide an essentially 5a signal for a Higgs-like resonance with 
mass of order 123-128 GeV. In the 77 final state, the ATLAS and CMS rates are roughly 1.9±0.4 and 1.6±0.4 times 
the Standard Model (SM) prediction. In the ZZ -t U channel, the ATLAS and CMS signals are roughly l.l±g;| and 
0.71q3 times the SM expectation, respectively. In the bb, t + t~ and WW — > Ivlv channels, the central value ATLAS 
rates are somewhat suppressed relative to the SM prediction but error bars are very large. The CMS signals in these 
latter channels are also somewhat suppressed and lie at least la below the SM prediction — no signal being observed 
in the t + t~ channel. Meanwhile, the CDF and DO experiments have announced new results [J] that support the 
~ 125 GeV Higgs signal and suggest an enhancement relative to the SM of the VF+Higgs with Higgs— > bb rate by a 
factor of 2 ± 0.6. 

Enhanced rates in the 77 channel have been shown to be difficult to achieve in the NMSSM [5], while remaining 
consistent with all relevant constraints, including those from LEP searches, _B-physics, the muon anomalous magnetic 
moment, a^ = (g — 2)^/2, and the relic density of dark matter, Vth 2 , when parameters are semi-unified at the GUT 
scale. By "semi-unified" we mean a model in which mo, mm, and Aq are universal at the GUT scale with NUHM 
relaxation for m 2 H , m\ and m s and general A\ and A K . Enhancements appear to be possible only if large values of 
the superpotential coupling A are employed and the a M constraint is greatly relaxed [B]. (See [7] for the first discussion 
of an enhanced 77 rate at large A in the NMSSM with parameters defined at the weak scale.) In this Letter, we 
pursue the case of generally large A and uncover a particularly interesting set of scenarios in which the two lightest 
CP-even Higgs bosons, hi and /12, both lie in the 123-128 GeV mass window. Phenomenological consequences are 
examined. 

For the numerical analysis, we use NMSSMTools 0(5] [ID] version 3.2.0, which has improved convergence of RGEs 
in the case of large Yukawa couplings and thus allows us to explore parameter regions that where left uncharted 
in [5]. The precise constraints imposed are the following. Our 'basic constraints' will be to require that an NMSSM 
parameter choice be such as to give a proper RGE solution, have no Landau pole, have a neutralino as the lightest 
SUSY particle (LSP) and obey Higgs and SUSY mass limits as implemented in NMSSMTools-3.2.0 (Higgs mass limits 
are from LEP, older TEVATRON, and early LHC data; SUSY mass limits are essentially from LEP.) 

Regarding B physics, the constraints considered are those on HR{B S — > X s j), AM S , AM d , BR(B S — > ^+/i~), 
BR(B+ t+v t ) and BR(5 X s ^t + ^) at 2a as encoded in NMSSMTools-3.2.0, except that we updated the 
bounds on rare B decays to 3.04 < BR(B S A s 7) x 10 4 < 4.06 and BR(£? ^ + /-0 < 4.5 x 10~ 9 ; theoretical 
uncertainties in B-physics observables are taken into account as implemented in NMSSMTools-3.2.0. 

Regarding dark matter constraints, we accept all points that have flh 2 < 0.136, thus allowing for scenarios in which 
the relic density arises at least in part from some other source. However, we single out points with 0.094 < ilh 2 < 0.136, 
which is the 'WMAP window' defined in NMSSMTools-3.2.0 after including theoretical and experimental systematic 
uncertainties. In addition, we impose bounds on the spin-independent LSP-proton scattering cross section implied 
by the neutralino-mass-dependent XenonlOO bound [IT] . (For points with ilh 2 < 0.094, we rescale these bounds by a 
factor of OAl/Qh 2 .) 

Our study focuses in particular on NMSSM parameter choices such that both and mh 2 lie within 123-128 GeV. 
We focus moreover on A > 0.1, a range for which it is known [B][7] that some enhancement, relative to the SM, of the 
Higgs signal in the 77 final state is possible. The degenerate situation is especially interesting in that an enhanced 
77 rate at ~ 125 GeV could arise as a result of the hi and /12 rates summing together, even if the individual rates 
are not full SM-like strength (or enhanced). 
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Above, we did not mention imposing a constraint on a M . Rough consistency with the measured value of requires 
that the extra NMSSM contribution, 5a M , falls into the window defined in NMSSMTools of 8.77 x 1CT 10 < Sa^ < 
4.61 x 10" 9 expanded to 5.77 x 1CT 10 < Sa^ < 4.91 x 1CT 9 after allowing for a la theoretical error in the NMSSM 
calculation of ±3 x 10~ 10 . In fact, given the previously defined constraints and focusing on A > 0.1, Sa^ is always too 
small, being at most ~ 2 x 10 -10 . Demanding 8a^ large enough to fall into the above window, or even come close to 
doing so, appears from our scans to date to only be possible if A < 0.1 5 a , for which the Higgs signal in the 77 and 
VV* (V = W, Z) final states for Higgs in the 123-128 GeV window is very SM-like. 

The main production/decay channels relevant for current LHC data are gluon-gluon and WW fusion to Higgs with 
Higgs decay to 77 or ZZ* — » 4£. The LHC is also beginning to probe W, Z+Higgs with Higgs decay to bb, a channel 
for which Tevatron data is relevant, and WW — > Higgs with Higgs— s- t + t~. For the cases studied, where there are 
two nearly degenerate Higgs bosons, we will combine their signals as follows in defining the mass and signal for the 
effective Higgs, h. First, for the individual Higgs we compute the ratio of the gg or VFVF-fusion (VBF) induced Higgs 
cross section times the Higgs branching ratio to a given final state, X, relative to the corresponding value for the SM 
Higgs boson: 

ff h, m _ iXffg -> hj) ERjhj -> X) 

T(gg -> h SM ) BR(h SM -> X)' { ' 

„ hi r( WW -> hj) BRjht ->■ X) 

VBF \ ) - T(ww ^ BR(hsM ^ x) , W 

where hi is the i th NMSSM scalar Higgs, and hsu is the SM Higgs boson. Note that the corresponding ratio for 
V* Vhi (V — W, Z) with hi —> X is equal to Ry BF (X). These ratios are computed in a self-consistent manner (that 
is, treating radiative corrections for the SM Higgs boson in the same manner as for the NMSSM Higgs bosons) using an 
appropriate additional routine for the SM Higgs added to the NMHDECAY component of the NMSSMTools package. 
Next, we compute the effective Higgs mass in given production and final decay channels Y and X, respectively, as 

m Y (X) = R y W TO/l i + R y ( x ) m h2 (gS 
~ R$(X)+R%?(X) 

and define the net signal to simply be 

R^(X) = R^(X) + R^{X). (4) 

Of course, the extent to which it is appropriate to combine the rates from the hi and /12 depends upon the degree of 
degeneracy and the experimental resolution. For the latter, we assume a TCS ~ 1.5 GeV H2]P] It should be noted that 
the widths of the hi and /12 are of the same order of magnitude as the width of a 125 GeV SM Higgs boson, i.e. they 
are very much smaller than this resolution. 

We perform scans covering the following parameter ranges, which correspond to an expanded version of those 
considered in [5]: < m < 3000; 100 < m 1/2 < 3000; 1 < tan/3 < 40; -6000 < A < 6000; 0.1 < A < 0.7; 
0.05 < k < 0.5; -1000 < A x < 1000; -1000 < A K < 1000; 100 < (i eS < 500. In the figures shown in the following, 
we only display points which pass the basic constraints, satisfy _B-physics constraints, have flh 2 < 0.136, obey the 
XENON100 limit on the LSP scattering cross-section off protons and have both hi and h^ in the desired mass range: 
123 GeV < m hl ,m, l2 < 128 GeV. 

In Fig. [TJ we display R^ (77) versus R^ (77) with points color coded according to m/, 2 — ■ The circular points 
have fi/i 2 < 0.094, while diamond points have 0.094 < flh 2 < 0.136 (within the WMAP window). We observe a large 
number of points for which m^, m/j 2 € [123, 128] GeV and many are such that ii^ (77) + (77) > 1. A few such 
points have Vth 2 in the WMAP window. These points are such that either ^^(77) > 2 or Rggdl) > 2, with the R 
for the other Higgs being small. However, the majority of the points with ii^ (77) + ^=(77) > 1 have Slh 2 below 
the WMAP window and for many the 77 signal is shared between the h\ and the hi- 

Based on these results, we will now combine the hi and /i2 signals as described above and present plots coded 
according to the following legend. First, we note that circular (diamond) points have £lh 2 < 0.094 (0.094 < f2/i 2 < 
0.136). We then color the points according to: 

• red for m/, 2 — < 1 GeV; 



The values for cr rcs quoted in this paper range from 1.39-1.84 GeV to 2.76-3.19 GeV, the better resolutions being for the case where 
both photons are in the barrel and the worse resolutions for when one or both photons are in the endcap. We anticipate that the more 
recent analyses have achieved substantially better mass resolutions, but details are not yet available. 
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FIG. 1. Correlation of gg — > (hi, hi) — > 77 signal strengths when both hi and /12 lie in the 123-128 GeV mass range. The 
circular points have flh 2 < 0.094, while diamond points have 0.094 < fl/1 2 < 0.136. Points are color coded according to 
mh 2 — rnh 1 as indicated on the figure. 



• blue for 1 GeV < ra^ 2 — mf ll < 2 GeV; 

• green for 2 GeV < m/ l2 — nif ll < 3 GeV. 

For current statistics and cr rcs > 1.5 GeV we estimate that the hi and hi signals will not be seen separately for 
m h2 - m h_,< 2 GeV. 

In Fig. 2 we show results for R h gg {X) with m h e [123, 128] GeV as a function of for X — 77, VV, bb. Enhanced 
77 and VV rates from gluon fusion are very common. The bottom-right plot shows that enhancement in the Wh 
with h — > bb rate is also natural, though not as large as the best fit value suggested by the new Tevatron analysis [4] . 
Diamond points (i.e. those in the WMAP window) are rare, but typically show enhanced rates. 

In Fig. pi we display in the left-hand plot the strong correlation between Rg g ("f^/) and Rg g (VV). Note that if 
Rg g (yf) ~ 1.5, as suggested by current experimental results, then in this model Rg g (VV) > 1.2. The right-hand plot 
shows the (anti) correlation between R gg (j"f) and Ry^,^ wh {bb) — Ry BF (bb). In general, the larger R gg {jj) is, the 
smaller the value of R^y»^ wh (bb). However, this latter plot shows that there are parameter choices for which both 
the 77 rate at the LHC and the W* — > Wh(—> bb) rate at the Tevatron (and LHC) can be enhanced relative to the 
SM as a result of there being contributions to these rates from both the hi and h 2 . It is often the case that one of the 
hi or hi dominates R gg (jj) while the other dominates R^-^whibb). This is typical of the diamond WMAP-window 
points. However, a significant number of the circular Qh 2 < 0.094 points are such that either the 77 or the bb signal 
receives substantial contributions from both the hi and the /12 (as seen, for example, in Fig. [I] for the 77 final state) 
while the other final state is dominated by just one of the two Higgses. We did not find points where the 77 and bb 
final states both receive substantial contributions from both the hi and hi- 

As noted above, there is a strong correlation between R^ g (^) and R gg (VV) described approximately by Rg 0(77) ~ 
1.25 Rg (VV). Thus, it is not surprising that the values for the gluon fusion induced 77 and VV cases are also 
strongly correlated — in fact, they differ by no more than a fraction of a GeV and are most often much closer, see 
the left plot of Fig. |4j The right plot of Fig. [4] illustrates the mechanism behind enhanced rates, namely that large 
net 77 branching ratio is achieved by reducing the average total width by reducing the average bb coupling strength. 

The dependence of R gg (jj) on A, n, tan/3 and fi e s is illustrated in Fig. [5] We observe that the largest Rggijj) 
values arise at large A, moderate k, small tan/3 < 5 (but note that R^ gg {a) > 1-5 is possible even for tan/3 = 15) and 
small /Lt e ff < 150 GeV. 
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FIG. 2. R h gg (X) for X = 77,VV,6&, and i? VBK (&&) versus mi,. For application to the Tevatron, note that Ry BF (bb) = 
Rw*^fWh{bb)- The color code here and in the following figures is green for points with 2 GeV < mi 2 — < 3 GeV, blue for 
1 GeV < m,h 2 — m-h! < 2 GeV, and red for rrih 2 — m; n < 1 GeV. 



Such low values of // e ff are very favorable in point of view of finetunig, in particular if stops are also light. Indeed 
a good fraction of our points with degenerate h\,h2 and #(77) > 1 features light stops with -Msusy = \J m lx m t 2 ~ 
1 TeV. The stop mixing is typically large in these cases, (A t — fi e fi cot/3)/MsuSY ~ 1-5-2. Moreover, the few points 
which we found in the WMAP window always have < 700 GeV. 

Implications of the enhanced 77 rate scenarios for other observables are also quite interesting. First, let us observe 
from Fig. [6] that these scenarios have squark and gluino masses that are above about 1.25 TeV ranging up to as high 
as 6 TeV (where our scanning more or less ended). The WMAP-window points with large Rg g (jj) are located at low 
masses of nig ~ 1.3 TeV and rriq ~ 1.6 TeV. 

The value of R h gg {ll) as a function of the masses of the other Higgs bosons is illustrated in Fig. [7J We see that values 
above 1.7 are associated with masses for the 02, /13 and H ± of order < 500 GeV and for the ai of order < 150 GeV. 
(Note that m Q2 ~ mh 3 — rn H ± ) While modest in size, detectability of these states at such masses requires further 
study. One interesting point is that m ai ~ 125 GeV is common for points with i?^ 9 ( 77 ) > 1 points. We have checked 
that #£1(77) is quite small for such points — typically < 0.01. 

Let us now focus on properties of the LSP. In the plots of Fig. [8j we display flh 2 and the spin-independent cross 
section for LSP scattering on protons, asi, for the points plotted in previous figures. We first note the rather limited 
range of LSP masses consistent with the WMAP window, roughly m^o g [60, 80] GeV. The corresponding ctsi values 
show a broader range from few x 10 -9 pb to as low as few x 10 -11 pb. Owing to the small /i e g, the LSP is dominantly 
higgsino, which is also the reason for Vth 2 typically being too low. The points with tth 2 within the WMAP window 
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FIG. 3. Left: correlation between the gluon fusion induced 77 and VV rates relative to the SM. Right: correlation between 
the gluon fusion induced 77 rate and the WW fusion induced bb rates relative to the SM; the relative rate for W* — > Wh with 
h — > bb (relevant for the Tevatron) is equal to the latter. 




FIG. 4. Left: effective Higgs masses obtained from different channels: mf l g ('yy) versus m^ 3 (VV). Right: 77 signal strength 
i4(77) versus effective coupling to bb quarks (C^) 2 . Here, C& 2 = [^(77)^' + i& 2 (77)<^ 2 ] / [^(77) + ^23(77)] • 



are mixed higgsino-singlino, with a singling component of the order of 20%, see the bottom-row plots of Fig. [8] 

It is interesting to note a few points regarding the GUT-scale parameters associated with the points plotted in 
previous figures. For the WMAP-window diamond points, m G [0.9,1.3] TeV, my 2 G [500,700] GeV, A G 
[-1.8,-1.0] TeV, A K G [-400,-250] GeV, A x G [-600,-400] GeV, m s (GUT) G [1.4,2.2] TeV, ms„ (GUT) e 
[2,2.2] TeV and m^f d (GUT) G [0.7,1.2] TeV; and, as shown in earlier figures, these diamond points have A G 
[0.58,0.65], k G [0.28,0.35], and tan/3 G [2.5,3.5]. Points with Rg g (ll) > 1.3 have m G [0.65,3] TeV, m 1/2 G 
[0.5,3] TeV, A G [-4.2,-0.8] TeV, A K G [-500, +450] GeV, A x G [-750, +550] GeV, m s (GUT) G [1.2,4.2] TeV, 
?7i H „(GUT) G [1.7,17] TeV, TO ff(J (GUT) G [-0,4.2] TeV, AG [0.33,0.67], k£ [0.22,0.36], and tan /? G [2,14]. 

We have already noted that it is not possible to find scenarios of this degenerate/enhanced type while predicting a 
value of <5a M consistent with that needed to explain the current discrepancy. In particular, the very largest value of 5a^ 
achieved is of order 1.8 x 10~ 10 and, further, the WMAP-window points with large R^ jg (^^f 1 VV) have da M < 6 x 10~ n . 

To summarize, we have identified a set of interesting NMSSM scenarios in which the two lightest CP-even Higgs 
bosons are closely degenerate and lie in the 123-128 GeV mass window. Large rates (relative to gg —> hsu 77 
or gg — > Hsm — > ZZ* — > At) for gg — > h\p —¥ 77 and gg — > hy2 ZZ* — > Al are possible, sometimes because 
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FIG. 5. Dependence of Rg g (-yy) on A, k, tan/3 and /x e ff- 



one of the rates is large but also sometimes because the rates are comparable and their sum is large. This suggests 
that, especially if enhanced rates continue to be observed in these channels, it will be important for the experimental 
community to be on the lookout for mass peaks in m 77 and UI41 that are broader than expected purely on the basis 
of the experimental mass resolution. In addition, the apparent mass in the 77 final state might differ slightly from 
the apparent mass in the A£ final state. Significant statistics will be required to resolve such features. 
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FIG. 8. Top row: flh and spin-independent cross section on protons versus LSP mass for the points plotted in previous figures. 
Bottom row: Q.h 2 versus LSP higgsino (left) and singlino (right) components. 



